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Abstract 



The understanding of oxygen reduction reaction (ORR) activity on perovskite oxide sur- 
faces is essential for promising future fuel cell applications. We report a comparative study of 
ORR mechanisms on LaBOj (S=Mn, Fe, Cr) surfaces by first-principles calculations based on 
density functional theory (DFT). Results obtained from varied DFT methods such as general- 
ized gradient approximation (GGA), GGA+U and the hybrid Hartree-Fock density functional 
method are reported for comparative purposes. We find that the results calculated from hybrid- 
functional method suggest that the order of ORR activity is LaMnOa > LaCrOa > LaFeOa, 
which is in better agreement with recent experimental results (Suntivich et ai, Nature Chem- 
istty 3, 546 (2011)) than those using the GGA or GGA+U method. 
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INTRODUCTION 



Oxygen reduction reaction (ORR), one of electrochemical energy conversion processes, plays a 
very important role in renewable energy technologies, including fuel cells and metal-air batteries. 
Searching for a highly active catalyst to replace noble metal cathodes is driven by the need of 
efficient and low-cost ORR catalysts, which are essential for mass marketing a fuel cell technology 
to address the world's energy needs. 

Recent experiments by Suntivich et al.—'^ shows that perovskite transition-metal oxides can 
exhibit high electrocatalytic activity for ORR in alkaline electrolytes. It is also suggested that 
ORR activity on perovskite oxide surfaces is related to the eg occupation in the B-site cation, 
which is indicative of the strength of bonding between transition metal ion and adsorbed oxygen. 
A moderate amount of eg-fiUing in perovskite oxides such as LaMnOs (eg = 1) yields higher 
activity as compared to other oxides with either too little (LaCrOs with eg = 0) or too much eg 
electron filling (LaFeOs with eg = 2). 

ORR on LaBOs perovskite oxide surfaces has been theoretically investigated extensively in 
the past using first-principles methods based on density functional theory.-"- All of these studies 
considered atomic and molecular oxygen adsorption on the LafiOs surface only for its application 
in solid oxide fuel cells. However, in alkaline fuel cells the ORR involves a more complicated 
reaction pathway, in which ORR intermediates such as hydroxides and peroxides- will form on 
the surface. Furthermore, in most studies the standard (semi) local DFT methods, local density 
approximation or generalized gradient approximation (GGA), are used, while only a few of them 
use GGA+U approach to describe strongly correlated electrons.-^ The role of GGA-i-i7 is to ad- 
dress on-site Coulomb interactions in the localized orbitals (such as d orbitals in transition metals) 
with an additional Hubbard-type term i7.— »^ It has been found that LaBOj surface energetics are 
strongly depend on the parameter U used in the calculation.- 

These recent studies further motivate us to use higher-level theoretical treatments, the so-called 
"parameter-free" hybrid functional approach, for a complete description of the ORR on the LafiOs 
perovskite oxide surfaces. A hybrid functional, a combination of exact nonlocal orbital-dependent 
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Hartree-Fock exchange and a standard local exchange-correlation functional, provides a signifi- 
cant improvement over the GGA description and enables accurate computation of electronic prop- 
erties and energetics of molecular systems as well as extended systems including transition metal 
oxides^"— without the need for system-dependent adjustable parameters or decisions of which 
electrons to localize. So far ORR activity on perovskite oxide surfaces is unexplored by hybrid- 
functional approaches. In the present work we aim at filling this gap. 

In this paper we report a comparative study of ORR activity on LaBOj (B=Mn, Fe, Cr) per- 
ovskite oxide surfaces using first-principles calculations. Calculations based on various DFT func- 
tionals such as GGA, GGA+U and hybrid functional are performed for comparative purposes. 
From binding energies of ORR intermediates, we obtain free energy changes at each step of the 
ORR and present overall free energy diagrams for each surface. The kinetics of ORR on each 
perovskite oxide surface are found to differ significantly with different DFT methods used in the 
calculation. We find that the hybrid functional method yields better agreement with recent experi- 
ments by Suntivich et al.^ while the results from the GGA and GGA-\-U methods fail to explain 
the experimental observations. 

Computational Methods 
DFT calculations 

Our calculations are performed using the plane-wave-basis-set Vienna ab initio simulation pack- 
age VASP— (version 5.2). The projector- augmented-wave (R^W) methods are used to describe 
the interactions between atomic cores and valence electrons, with a kinetic energy cutoff of 500 eV 
employed in all simulations. The 3p-semi-core states are treated as valence states for transition 
metal atoms Cr {?ip^?id^As^), Mn (ip^'id^As^), and Fe (3p^3J^45^). For the exchange-correlation 
functional we used the GGA method with Perdew-Burke-Ernzerhof (PBE) formulation. We also 
apply the GGh+U method to reduce the self- interaction error and improve the description of corre- 
lation effects. For the GGA+U calculation we make use of the standard Dudarev implementation— 
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where the on-site Coulomb interaction for the localized orbitals is parametrized by Ueff = U — J 
using the PBE functional. We apply the optimized effective interaction parameter Ueff for the 
metal atoms in LaBOs (Ueff = 4, 4 and 3.5 eV for Mn, Fe and Cr, respectively), determined by 
fitting the enthalpies of the oxidation reactions. These values have previously been shown to 
provide a description of LafiOs electronic structure that is in good agreement with the available 
experimental data.- Furthermore, we perform single-point total energy calculations with the hybrid 
functional approach developed by Heyd-Scuseria-Ernzerhof (HSE06)— for the approximation of 
the exchange-correlation energy and potential. In the HSE06 approach, one quarter of the PBE 
short-range exchange is replaced by the exact Hartree-Fock exchange, and the full PBE correlation 
energy is included. The range-separation parameter is set to be 0.2 A \ 
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Figure 1: (Color online). Optimized geometry of bare LaMnOs (001) surface (a) and surfaces with 
adsorbed ORR intermediates 00* (b), O* (c), HO* (d), and HOO* (e) intermediates. The dashed 
lines are the boundaries of the supercell. 

The LafiOs perovskite oxide surfaces are simulated with a slab including nine atomic layers 
with transition metal 5-terminated (001) surfaces at each side. Such surfaces are chosen simply be- 
cause the (001) surfaces are generally the most stable in perovskites and the redox active transition 
metal atoms are expected to be responsible for ORR catalytic activity.- The distorted perovskite 



5 



structure of pnma is considered to be the structure identified in the experiment of Ref . (author?) ^ , 
and the equilibrium bulk geometry with optimized lattice constants is used to construct the LaBOs 
slab geometries. A large vacuum spacing of 15 A was used between two slabs to prevent interaction 
between the system with its images in adjacent unit cells. We use the 1x1 surface unit cell, and 
two-dimensional periodic boundary conditions are applied for the surface directions with a 6 x 6 
Monkhorst-Pack type of fc-point sampling. The four bottom atomic layers of the slab are fixed at 
the optimized bulk lattice constant, while the top five layers as well as the ORR intermediate on the 
surface are fully relaxed. In all calculations a dipole correction is applied. The ORR intermediates 
are modeled by the same slab with the intermediate adsorbed at the B sites on the top surface of 
the slab. Adsorbates placed above the surfaces correspond to a coverage of 50%. The geometries 
are optimized until the force on each atom falls below the convergence criterion of 0.02 eV/A. The 
optimized structures of bare LaBOs surfaces and surfaces with adsorbed ORR intermediates 00*, 
O*, HO* and HOO* (asterisk denoting adsorbed species hereafter), taking LaMnOs as an exam- 
ple, are shown in Fig. \\\ (a)-(e). For the HSE06 calculations, four bottom atomic layers of each 
La503 slab are removed in order to save computational effort. 

The magnetic structures of LafiOs perovskites are complicated due to phase transitions at dif- 
ferent temperatures. For simplicity and consistency we apply the ferromagnetic ordering for all 
the studied LaBOs perovskites, by an initial assignment of non-zero (parallel) spin to all the B 
atoms. The self-consistent procedure determining the electron and spin density distribution for 
each geometry are always carried out leaving the total spin unconstrained and free to evolve. We 
expect our assumption will not introduce large errors in our energetic results, as we focus on the 
adsorption energy and relative free energy of the ORR intermediates on the surface, in which the 
total energy difference caused by changing the relative spin orientations inside the slab should have 
negligible effect. 
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Free energy diagrams for oxygen reduction reaction 

The following four-electron ORR reaction pathway for LaBOs at the cathode in alkaline elec- 
trolytes has been suggested- (overall process O2 + 2H2O + 4e^ — )• 40H^): 

1) surface hydroxide displacement as: 

La503-OH + O2 + e — > LafiOa-OO + OH" (1) 

2) surface peroxide formation as: 

La503-00 + H2O + e — > LaBOs-OOH + OH (2) 

3) surface oxide formation as: 

La503-OOH + e~ — > LafiOs-O + OH" (3) 

and 4) surface hydroxide regeneration as: 

La503-0 + H2O + — > LafiOa-OH + OH (4) 

The free energy change of each ORR reaction step is calculated based on a computational 
hydrogen electrode model suggested by N0rskov et al?^ This method has been shown to predict 
trends for the oxygen reduction reaction on metals quite well.— In this method the potential affects 
the relative free energy through the chemical potential of the electrons in the electrode. We convert 
the calculated DFT energies into Gibbs free energies by adding entropic {TS) and zero-point energy 
(ZPE) corrections to the ORR intermediates, so that 

AG = A£ + AZPE-rA5 + AG<i> (5) 
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where is the calculated DFT reaction energy, AZPE is the change in ZPE and AS is the change 
in the entropy. ZPE corrections and entropies of the ORR intermediates are calculated from the 
vibrational frequencies according to standard methods, and those of the gas-phase molecules are 
obtained from thermodynamics databases. AG^ is the effect of electrode potential which is ap- 
plied by shifting the free energy change AG by AG^ = e4>, where e is the elementary charge 
and 4> is the potential difference between electrode and counter electrode (versus reversible hy- 
drogen electrode, RHE). The equilibrium potential 4>eq corresponds to zero net reaction free en- 
ergy (E/Li Gi = 2AGw + 4e4> = 0) of the overall ORR process, thus we have 4>eq = -AGw/(2e), 
where AGw = G(H20) — G(H2) — G(02) /2 is the free energy of formation of water from H2 and 
5O2. The ORR potential 4>orr corresponds to the highest potential 4> at which all steps along 
the reaction decrease the free energy. The theoretical ORR overpotential is then calculated by 

77 =4>eq-4>ORR. 

Furthermore, effect of water on the ORR activity is also taken into account by a solvation 
correction to ORR intermediates, since reactions occur in the presence of water in electrolytes .- 
The solvation corrections are obtained from a previous study of platinum-catalyzed ORR using 
the Poisson-Boltzmann implicit continuum model. The solvation correction energies are —0.32, 
-0.47, -0.75 and -0.54eV for 00*, O*, HO*, and HOO* intermediates, respectively. 

Results and discussion 

Surface Energetics and Electronic properties 

We first calculate the DFT binding energy for each ORR intermediate on the LaBOs (001) surfaces. 
The DFT binding energies of the ORR intermediates are calculated according to A^adsorbate = 

£'LaB03-adsoibate — £^LaB03 — ^'adsorbate » where ^'LaBOs-adsorbate; ^'LaBOj and ^adsorbate are the energies 

of the surface with adsorbed intermediate, the bare surface and the isolated ORR intermediate, re- 
spectively. The negative sign of ^adsorbate corresponds to energy gain of the system due to adsorp- 
tion of ORR intermediate, and the more negative the ^adsoi-bate value the stronger is the chemical 
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Table 1: Calculated DFT binding energy of each ORR intermediate on LaBOs surfaces 
(unit: eV). The DFT binding energies of the ORR intermediates are calculated according to 

^^adsorbate ~ ■£^LaB03-adsorbate - ^LaBOj - ^adsorbate. The absorbtion energies AE^^^^^^^^^ relative 
to energies of H2O and H2 are also given (in parentheses, see text for definition). 



siruciurca 


IIICIIIUU 




A 


AF^^TT fAF^^ \ 
i^^OOn V^^OOH-' 






GGA 


-2.927 (0.417) 


-0.659 


-1.424 (3.739) 


-3.755 (1.913) 


LaMn03 


GGA+U 


-2.107 (1.236) 


-0.205 


-0.724 (4.439) 


-1.987 (3.681) 




HSE06 


-1.647 (1.355) 


0.162 


-0.482 (4.844) 


-1.278 (4.015) 




GGA 


-2.686 (0.658) 


-0.491 


-1.224 (3.939) 


-3.457 (2.211) 


LaFeOs 


GGA+C/ 


-2.287 (1.057) 


-0.199 


-0.794 (4.369) 


-2.134 (3.534) 




HSE06 


-1.778 (1.223) 


0.328 


-0.028 (5.298) 


-1.368 (3.926) 




GGA 


-3.276 (0.068) 


-0.931 


-1.786 (3.377) 


-4.986 (0.682) 


LaCrOs 


GGA+t/ 


-2.429 (0.915) 


-0.172 


-1.014 (4.149) 


-3.515(2.153) 




HSE06 


-1.858 (1.144) 


0.305 


-0.678 (4.648) 


-2.972 (2.322) 



interaction between the adsorbate and the surface. It should be also noted that our calculations give 
the reference energies of —9.86 and — 17.04eV for O2 with the GGA and the hybrid functional, 
respectively, and no correction for this has been made in the values of this work. While such a 
correction would alter the absolute binding energies it has no effect on the relative difference in 
binding energies between surfaces. The results are summarized in Table[I]for comparison among 
different DFT functional. It is clear that for each ORR intermediate the binding energy (abso- 
lute values) decreases from the GGA value to that of GGA+U and then especially strongly to the 
HSE06 hybrid-functional result. This holds for all three perovskite surfaces. 

Among the four ORR intermediates, the O* binding energy has the strongest dependence of the 
calculating method. Taking LaMnOs as an example, the GGA binding energy Eq is about 1.8eV 
larger (in its absolute value) than that of GGA-i-i7, and it is about 2.5 eV larger than the HSE06 
value. For the binding energy £00 ^ relatively weaker dependence can be found. For LaMnOs 
the difference between calculated values of £00 from the GGA and GGA+U methods is about 
0.45 eV, and that between the GGA and HSE06 methods is about 0.8eV. Similar behaviors can 
be also found for LaFeOs and LaCrOs in Table \T\ These can be explained by the fact that the 
interaction between atomic oxygen and the surface transition metal B ion is much stronger than 
the other ORR intermediates, thus the effect of correlation corrections in the the GGA-i-i7 and 
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HSE06 calculations are much more obvious on the O* binding energies. The calculated charge 
transfer results further support this point. In Table [2l we show the effective Bader charges (based 
on the real-space-charge density of transition metal ion and adsorbed ORR intermediate on 
each LaBOs surface. It is clear that the charge transfer between the surface and the adsorbed ORR 
intermediate is large for LafiOs-O but much smaller for LafiOs-OO. Our results are consistent with 
a recent ab initio GGA+U study by Lee et al. , ^ which also shows that the O* binding on the LaBOs 
surface has a stronger U^ff dependence than the 00* binding. In our calculations, the even larger 
correction of binding energies in the HSE06 hybrid functional results compared to the corrections 
in the results from the GGA+U method is primarily due to the fact that the f/gff is applied only at 
the transition metal B atoms in the GGA+U calculations. 

Table 2: Effective Bader charges q (in e) of transition metal ions and adsorbed ORR inter- 
mediates (in parentheses) on LafiOs surfaces. 



perovskites 


method 


bare surface 


La503-00 


La503-OOH 


La503-0 


La503-OH 




GGA 


1.66 


1.78 (-0.31) 


1.81 (-0.36) 


1.80 (-0.51) 


1.82 (-0.39) 


LaMnOs 


GGA+U 


1.71 


1.76 (-0.16) 


1.81 (-0.30) 


1.86 (-0.52) 


1.89 (-0.44) 




HSE06 


1.83 


1.91 (-0.25) 


1.96 (-0.37) 


1.93 (-0.46) 


1.99 (-0.42) 




GGA 


1.57 


1.64 (-0.24) 


1.68 (-0.33) 


1.70 (-0.52) 


1.74 (-0.47) 


LaFeOs 


GGA+U 


1.79 


1.84 (-0.15) 


1.85 (-0.35) 


1.73 (-0.54) 


1.83 (-0.47) 




HSE06 


1.81 


1.94 (-0.18) 


1.93 (-0.41) 


1.87 (-0.57) 


1.95 (-0.55) 




GGA 


1.78 


1.92 (-0.40) 


1.91 (-0.39) 


2.01 (-0.55) 


1.93 (-0.42) 


LaCrOs 


GGA+U 


1.81 


1.97 (-0.18) 


1.93 (-0.29) 


2.02 (-0.51) 


1.99 (-0.39) 




HSE06 


1.94 


2.02 (-0.28) 


1.91 (-0.39) 


2.17 (-0.56) 


2.14 (-0.43) 



In comparing the calculated binding energies of the three LaB03 systems, we notice that for 
each ORR intermediate the binding energy (in its absolute value) is generally larger for LaCr03 
than that of LaMn03 or LaFe03, with the only exception the 00* intermediate with the GGA+U 
and HSE06 results. In particular, the binding energy Eq for LaCrOs is (depending on the method) 
about 1.2~1.7 and 1.4~1.6 eV more negative than that of LaMnOs and LaFe03, respectively, 
indicating a stronger interaction between the LaCr03 surface and the atomic oxygen adsorbate. 
For the 00* intermediate, we find that though the GGA results still show a ~ 0.3 eV larger binding 
energy (in its absolute value) for LaCr03 as compared to the LaMnOs (~ 0.4eV larger as compared 
to LaFeOs), the GGA-i-i7 method yields very similar binding energies of molecular oxygen for all 
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three perovskite surfaces. The HSE06 method, however, gives the largest binding energy £00 with 
a positive sign for LaFeOs surface, indicating a relatively weak binding ability of LaFeOs surface 
for 00* intermediate. This is also true for the HOO* intermediate where the HSE06 method yields 
smaller binding energy ^ooH (in its absolute value) for LaFeOs than the other two perovskites. 

Li Table [U we also show the calculated adsorption energies of HO*, HOO* and O* intermedi- 
ates relative to H2O and H2, according to the following reactions as defined in Ref. (author?)—: 

^OH = £^LaB03-0H " ^'LaBOs " (^HiO " jE^l^) (6) 
^OOH = £'LaB03-00H -£'LaB03 - (2£'h20- l^'Ha) (7) 
A£'o'^ = ^LaBOa-O — £'LaB03 — (^HaO — ^Hi) • (8) 

where En^o, and En^ are the calculated energies of H2O and H2 molecules in the gas phase, 
respectively. Our calculated adsorption energies of HO*, HOO* and O* intermediates with the 
GGA-PBE method are very close to the recent reported results by Man et ah— using the GGA- 
RPBE method, with less than 0.1 eV difference for LaFeOs and LaCrOs. For LaMnOs our GGA 
adsorption energies are slightly larger than the results reported in Ref. (author?) , with 0.55,0.15 
and 0.27 eV difference for the O*, HO* and HOO* intermediates, respectively. The discrepancy 
may be because of different pseudopotentials (ultrasoft) used in the Ref. (author?)^-. 

ORR on La503 surface 

The calculated free energy change of each ORR step at zero electrode potential (4> = V vs. RHE) 
are listed in Tabled AGi, AG2, AG3, and AG4 correspond to the free energy changes in the ORR 
reaction steps of Eqs. (1), (2), (3) and (4), respectively. All reaction steps are exothermic for 
all three perovskite oxide surfaces with different calculation methods. However, for each surface 
GGA+U and HSE06 methods give quite different relations between different reaction steps respect 
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Table 3: Calculated free energy change AG for each ORR step in LafiOs (unit: eV). * denotes 
the potential-determining step for the ORR with the smallest |AG|. 



structures 


method 


AGi 


AG2 


AG3 


AG4 




GGA 


-0.642* 


-0.649* 


-2.370 


-1.000 


LaMnOs 


GGA+t/ 


-1.008 


-0.403* 


-1.302 


-1.948 




HSE06 


-0.758 


-0.526* 


-1.372 


-2.165 




GGA 


-0.715 


-0.618* 


-2.272 


-1.057 


LaFeOs 


GGA+U 


-0.822 


-0.479* 


-1.379 


-1.981 




HSE06 


-0.462 


-0.237* 


-1.916 


-2.206 




GGA 


-0.564 


-0.739 


-3.239 


-0.119* 


LaCrOs 


GGA+t/ 


-0.652* 


-0.726 


-2.540 


-0.743 




HSE06 


-0.405* 


-0.864 


-2.870 


-0.682 



to the GGA values. The size of the ORR potential-determining step can be estimated according to 



G = min(|AGi|, IAG2I, IAG3I, IAG4I). (9) 



In the following we further analyze the catalytic performance by calculating the free energy dia- 
grams under equilibrium and ORR potentials for each LaBOs surface. 

Fig. [2I shows the free energy diagrams of the ORR on the LaMnOs (001) surface with different 
calculation methods for representative potentials 4>eq and 4>orr. At the equilibrium potential 4>eq, 
for the GGA method it is clear that only reaction step (3), surface O* formation, are energetically 
downhill (exothermic) and all other three steps are uphill (endothermic). The O* intermediate is 
stabilized on the LaMnOs surface. This is very similar to the case of ORR on Pt(lll) surfaces 
studied by N0rskov et al.^-'^^ The first two steps of the ORR process, the surface 00*/H0* ex- 
change and the HOO* formation, are the two potential-determining steps. The fact that these two 
steps become exothermic at nearly the same ORR potential indicates that LaMnOs has a reactivity 
close to optimal, as the free energy changes in the four reaction steps are linearly related; thus 
increasing one of the step will decrease another. The theoretical ORR overpotential is 0.52 V, 
which is very close to the theoretical overpotential of 0.48 V for ORR on Pt(l 1 1).^^ In contrast to 
the strongly bound O* intermediate with the GGA functional, applying the GGA+U and HSE06 
functionals change the results by giving much higher relative free energy of O* intermediate in the 
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HO* 00* HOO* O* HO* 

Reaction Coordinate 



Figure 2: (Color online). Free energy diagrams of the ORR on the LaMnOs (001) surface. Solid 
and dashed lines represent reactions at the ORR potential (4> = 4>orr) and the equilibrium potential 
(4> = ^eq)^ respectively. The ORR potential corresponds to the highest potential 4> at which all 
steps along the reaction decrease the free energy. The theoretical overpotentials r] are marked for 
results with different functionals. 

ORR process. This should mainly be attributed to the largely decreased interaction strength be- 
tween the O* intermediate and the LaMnOs surface for the GGA+U and HSE06 functionals with 
respect to the GGA result, consistent with the calculated binding energy as shown in Tabled] The 
HO* intermediate becomes the most stable one on the LaMnOs surface and the surface hydroxide 
regeneration (step 4) becomes exothermic. The potential limiting steps lie in the second electron 
transfer step of the ORR (OOH* formation), and the theoretical ORR overpotentials are 0.76 V 
and 0.68 V for the GGA+U and HSE06 results, respectively. 

The calculated free energy diagrams for ORR on the LaFeOs (001) surface under 4>eq and 
^ORR are shown in Fig.|3l It has previously been suggested that for LaFeOs the surface 00*/H0* 
exchange does not gain sufficient energy, and thus the ORR kinetics are limited by the rate of 
00*/H0* exchange (step 1).^ However, our calculation indicates that this reaction step is only 
the secondary potential-determining step. The primary potential-determining step is located at the 
HOO* formation (step 2), which has a slightly larger free energy change than that of the 00*/H0* 
exchange step under equilibrium potential independent of the choice of the calculation method. 
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Figure 3: (Color online). Free energy diagrams of the ORR on the LaFe03 (001) surface. 



The GGA calculation for LaFeOs yields results very similar to the calculated free energy diagram 
for LaMnOs, for which we obtain a low ORR overpotential of 0.55 V. Also similar to LaMnOs, 
the self-interaction corrections in the GGA+U and HSE06 calculations lower the binding energy 
of O* on the LaFeOs surface and give clearly higher relative free energy for the O* intermediate 
than the HO* and 00* intermediates in the ORR process. Moreover, for the HSE06 method the 
free energy change in the two endothermic steps, 00*/H0* exchange and HOO* formation, are 
increased as compared to the results from GGA and GGA+U. As a result, while the GGA+U 
method gives a theoretical ORR overpotential of 0.69 V, comparable to the GGA method, the 
overpotential calculated for the HSE06 hybrid-functional is 0.97 V, which is almost twice as large 
as that of the GGA value. 

For ORR on a LaCrOs (001) surface, at the equilibrium potential the O* intermediate is 
strongly stabilized on the surface for the GGA functional, as shown in Fig. HI This situation 
does not change much even if we apply yhe GGA+U method or HSE06 functional in the cal- 
culation. For the GGA functional it is clearly shown that the surface HO* regeneration reaction 
(step 4) is the potential-determining step. The corresponding ORR overpotential is 1.05 V, much 
higher than that of LaMnOs or LaFeOs. This agrees very well with the experimental observations 
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Figure 4: (Color online). Free energy diagrams of the ORR on the LaCrOs (001) surface. 



in Ref. (author?)-, which suggests that the O* intermediate on the LaCrOs surface is not suf- 
ficiently destabilized and the ORR kinetics are limited by the rate of surface OH* regeneration. 
However, using both the GGA+U method and HSE06 functional we find the relative free energy 
of O* intermediate becomes about 0.6eV higher with respect to the GGA result, which makes the 
HO* regeneration step no longer the most endothermic step in the ORR. Instead, the H0*/00* 
exchange step (step 1) becomes the potential-determining step. For the HSE06 functional the 
calculated ORR overpotential is 0.80 V, which is 0.12 V larger than that of LaMnOa but 0.17 V 
smaller than that of LaFeOs. The overpotential of 0.51 V calculated from the GGA+U method is 
much lower than that of LaMnOs (0.76 V) or LaFeOs (0.69 V), in conflict with the experimental 
observations in Ref. (author?)-. We therefore conclude that the GGA-i-i7 method with U^ff applied 
on the Cr d state of LaCrOs may not be a physically appropriate approach to calculate the surface 
reaction energies. 
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Summary 

In summary, we conduct a comparative first-principles study of the ORR activity of the three 
perovskite oxides LaMnOs, LaFeOs and LaCrOs. In addition to the extensively used GGA and 
GGA+U methodologies, we also apply the advanced hybrid-functional method. We find that the 
calculated surface binding energies of ORR intermediates are strongly dependent on the method, 
and the free energy diagrams of ORR are described quite differently in GGA, GGA-i-i7 and hybrid 
functional approaches, especially for LaFe03 and LaCr03 surfaces. 

We show that Cr-sites on the LaCr03 surface are better adsorption centers for atomic oxygen 
than the other two transition metals, but the strong binding of the O* intermediate does not fa- 
vor the surface hydroxide (OH*) regeneration step of the ORR process. For the LaFeOs surface 
the interactions between the Fe-site and the 00* and HOO* intermediates are relatively weak, 
and the calculated ORR free energy diagram indicates that the HOO* formation and 00*/H0* 
displacement reactions are the primary and secondary potential-determining steps, respectively. 
These findings agree well with much of the gg-fiUing model deduced from experiments. The 
results calculated from the hybrid-functional method suggest that the order of ORR activity is 
LaMnOs > LaCrOa > LaFeOs. This is in better agreement with recent experimental observation 
in Ref . (author?) ^ than those from the GGA or GGA+U method. The GGA results yields similar 
free energy diagrams and ORR activities for LaMnOa and LaFeOa, while the GGA+U suggest that 
LaCrOs has the lowest ORR overpotential. In neither case does the obtained order of ORR activity 
agree with the experiment. 
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